iU3DXTXasr and StraTRZkCOrXON TEST CASES ^ 
AND HB7HOD 2^ SYST^H THEXtHFOR 

FXBZiD OF THE XKVENTXON 

The present Invention relates to the testing of 
floating-point arithmetic units, and^ more particularly, 
to the generation of numerical test cases for binary 
floaring-poini: adders- 

BSkCKSBOTmO OF THE XI3VEKTXON 

When developing integrated circuits that perform 
floating-point arithmetic^ designers typically base the 
representations of floating-point (FP) numbers and the 
constraints on the results of arithmeric operations on 
published standards, such as the well-^known ^^lEEE 
standard for binary floating point arithmetiCy An 
American National Standard, ANSI/IEEE 2td 754-1985^', 
Adherence to such standards guarantees that the circuitry 
will perform floating-point arithmetic with acceptable 
and predictable results. Although it is a relatively 
straightforward task to implement floating-point 
standards in a floating-point hardware unit, designers 
usually make modifications in the implementation to 
Improve performance in special cases. Because of this, it 
is necessary to xj^erify compliance of the finished design 
to the selected standard. In many instances, errors in 
floating-point implementation escape detection and find 
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rheir way into production. Cases such as the well— knoxrjn 
Pentium bug show that the verification process in this 
area is still far from being optimal. The ever-growing 
demand for increased performance, reduced time-to-market, 
and decreasing tolerance for errors all combine to make 
verification increasingly harder. The term 

'^floating-point unit'' herein denotes any device or system 
capable of performing binary floating-point computations 
by any means including, but not limited to, hardware^ 
firmware, programmable logic arrays / and software. 

There are many places problems can occur In the 
implementation of a floating-point unit, ranging from 
data problomg on single instructions to the correct 
handling of sequences of instructions in which 
back-to~back ©vents challenge superscalar 

implementations • This complexity stems both from the 
interpretation of the specification (architecture) as 
well as the peculiarities of the iinplementation 
(microarchitecture) . 

Although there is on-going work to develop formal 
proofs of adherence to a standard, formal methods are 
still far from providing a complete answer to the 
problem. The simulation of test cases (generating the 
test case data, running the test case on the target 
floating-point unit, and confirming the accuracy of the 
result) , has traditionally been employed for 
verification, and therefore remains the foundation of the 
verification process. 

It is the generating of floating-point test cases 
that is of interest regarding the present invention. 
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Yesr Cas^ Generarinq Background 

First, it: is clear that there is an enormous ^ 
practically unlimited number of different calculation 
-to test. In practice tlien^ simulation can be done 
on only a ve^ry small portion of the ©Kisting space. 
Reducing the enormous nuinber of potem:ial test cases to a 
manageable number that can actually be tested is done 
through placing suitable conszralnts on the machine 
numbers so that the constrained set o£ machine numbers 
used in the test will be representative of a particular 
aspect of testing, but will still constitute a sufficient 
nuinber of cases for thorough testing. (Constraints are 
discussed below,) 

The rationale behind verification-by-simulation is 
that one acquires confidence in th© correcnness of a 
floating-point unit design by running a set of test cases 
that encompass a sufficiently large number o£ different 
cases^ which in some sense is assumed to be a 
representative sample of rhe full space- The abiliry of 
the floating-point unit design to correctly handle all 
cases ±s inferretl from the correcr handling of the cases 
actually tested. 

To confidently make the above inference requires the 
building of a set of test cases that covers all special 
implementations of the floating-point unit design. The 
problem then becomes one of how best to do this. Since 
both the architecture specification and the 
microarchitecture implementation tend to yield a myriad 
of special cases^ generating the test cases using a 
uniform random distribution over the entire 
floating-point space would be highly inefficient. For 
IL 9-2001-0023 3 



example^ Ir is coimnon zhat executing an BM)D instruction 
that results in a sum of zero exercises a specific part 
of the design logic^ and therefore such a case should be 
verified. The probability of randomly generating two 
floating-point numbers that add to zero, however, is 
extremely low. Therefore^ prior-art random test 
generators usually possess some internal Testing 
Knowledge (TK) to bias the test generation towards cases 
of interest. Such test generators are described in 
''Model-Based Test Generation For Processor Design 
Verification'^ by Y.Lichtenstein, Y. Malka and A. Aharon, 
Innovative Applications or Artificial Xntelllg^nae 
(lAAI), A?^I Press, 1994; ^^Constraint Satisfaction for 
Test Program Generation" by Fournier, D. Lewin, 

Levinger, E. Roytman and Gil Shurek, Int. Phoenix. 
ConrejT&nce on Computers and Communications, March 1995,- 
and ^^Test Program Genera tion for Functional Verification 
of PowerPC Processors in IBM" by A. Aharon, Goodman, 
M. Levinger, Y, Lichtenstein, Y. Malka, C. Metzger, M. 
Molcho and G, Shurek, 32nd D&slgn Automation Conference, 
San Francisco, June 1995, pp. 279-285, 

In effect, TK changes the probability distribution 
of a test space, better adapting that test space to 
existing knowledge. In a test-generator described in the 
foregoing references^ the TK is in the form of C-language 
functions (called ^'generation functions'') which can be 
added incrementally to the generator, such as by the 
users themselves. A serious limitation of this prior art 
approach, however, is that such generation functions are 
very complex and difficult to write, requiring a deep 
understanding of the Floating Point unit design. In 
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practice, rhen, very few generating functions havB been 

added. 

Definitions 

The following rerms and abbreviations axe used 
5 herein: 

s±&s€^ esposenfe rhe siom of the eKponent and a 

constant "bias'''' selected so thar the 
biased exponent is always non-negative 
over the exponent's range. The use of 
a biased exponent allows representing 
both positive and negative exponents 
without requiring a sign* 
(a ^"denormalized nuiriber''') a non-zero 
FP number whose exponent has a 
reserved value (usually the mlnliuuxn 
pe?rmitted by the format) , and whose 
explicit or implicit leading bit of 
the significand is zero, 
the component of a binary 
floating-point number that normally 
specifies the integer power of 2 which 
is multiplied by the significand to 
esspress rhe value of the represented 
number. In certain formats, reserved 
values of the exponent are used to 
indicate that the number is 
denormalized. 

Floating-Point- A binary loating-point 
number contains a sign^ an exponent/ 
and a significand- 
Fia^ied-Point 
5 
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Ma,sk a teiaplate for a binary number. A mask 

is a series of characters r each of 
which represents allowable values for 
the bit in the corresponding position 
5 of the binary number- There is exactly 

one character in a mask for each bit 
o£ the corresponding binary number^ 
with characters in a one-to-one 
correspondence with the bits. 
Allowable mask characters include ^0', 
'1', and ^x' . A ^0' characrer 
indicates that the corresponding bit 
of the binary nTimber must be a 0^ 
whereas a ^1' chargicter indicates that 
^5 the corresponding bit of the binary 

number must be a 1. An character 
is a ^Mon't care", meaning that the 
corresponding bit of the binary number 
may be eirher 0- or 1. in a 
floating-point number, it is possible 
ro specify separate masks for sign, 
eKponent/ and significand. 
^«wr (""Not a Number''} a symbolic entity 

encoded in FP format. 
2S Noxmal (a ""normalized number") a non-zero FP 

number whose explicit or Implicit 
leading bit of the significand is one. 
3±gxx£± ^^n < S the component of a binary FP number 
that consists of a single explicit or 
inrplicit leading bir to the left of an 
implied binary point and a field of 
IL 9-2001-0023 6 
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fracrion bits ro the right of the 
binary point. Even in implementations 
where the signif icand' s leading bit 
(to the left of tbe binary point) is 
iit^lied and is not eKplicitly present^ 
the eignificanca still includes that 
leading bit. Note that a signifioand 
differs from a "mantissa^' in that a 
mantissa includes only the fraction 
field to the right of the binary 
point, whereas a signifioand also 
includes the bit to the left of the 
binary point. The concept of mantissa 
may be defined in such a way to 
include all of the bits of the 
signifioand, including the leftmost 
bit and not only th© fraction bits. In 
this case and for a binary normalized 
Signifioand, the mantissa would simply 
be the significand shifted one place 
to the right, or divided by a scaling 
factor of 2. 
The Set of Machi ne Murobers 

For purposes of illustration, a non-limiting enample 
of a binary floating-point number system is the IEEE 
standard 754 previously referenced. We assume that three 
integer constants are given, p. ^^^^^^ 

numbers are those which can be represented in the form 
(-1)" X 2^ X bobaba...bp-i, where ^ e {o, 1) represents 
rhe Sign or v. e, representing the exponent of v. is an 
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integer satisfying B^iti ^ -e: ^ JSiaa*, Ttie bit valuess are 
denoted as Jbi ^ (0, 1}, and p is the ^'precision" of the 
system. The significand is bobibz- * "i^xr whose binary 
point lies between Jbo and 2ji. All machine niimbers v that 
satisfy \v\ S are assumed to be normalized {Jbo D - 

Those machine numbers which are smaller in magnitude than 
2smdn (including zero) have E - S^pin and are denormalized 
(ho — 0) - Thus,r each macnine niimber has a unique 
representation (note that the IBEE standard 754 requires 
the same uniqueness for single and double formats but not 
for extended formats) • 

Binary Representations of Machine Numbers and the Mask 
Constraint 

Machine numbers are herein represented as ^strings of 
binary digits (bits) . This is true for fixed point 
numbers as well as for floating point numbers, A mstsk 
related to a number is a string of characters of the same 
length (number of bits) as the number^ all of whose 
characters are in the s^t {^0'^ } . A number and a 

mask are compatible if all the following conditions are 
met: 

• The number and the mask are of the same length 

(hence each bit of the number corresponds to a 
specific character of the mask) ; 

• For each '^l' character of the mask, there is a 1 
value in the corresponding bit of the number; and 

• Per each ^0' character of the masJc, there is a 0 
value in the corresponding bit of the number. 

If one or more of the above conditions are not met, 
the number and the mask are i/^aomp^tlble. Thus, a ^1' or 
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a ^0' character of the mask determines uniquely the value 
of the corresponding bit of the nufinber* An character 
in the mask leaves the corresponding bit value of the 
number under ermined- 

A number is constrained by requiring it ro be 
compatible wirh a given mask. 

Where it is not convenient to represent a 
f loaring-point machine number by a single string of bits^ 
it is possible {although not necessary) to split such a 
representation into a triplet of numbers: 

Sign: A string of one bit/ which is 0 for a 

* + ' and 1 for a . The numerical 

value is denoted by 5 (s « 0 or 1) . 
Biased exponent: A string of ?/ bits. This is 
interpreted to be a binary Integer with 
the numerical value where 0 < e :^ 

2*^-1. In a non-limiting generalisation 
from the single and double foannats of 
the IEEE standard 754, Eain ^ 2 - Z'^S 

-Eiaax = bias ^ 2""^ - 1. 

Signifieand: A string of p bits^ iSoi^xJba i^i- Unlike 

the single and double formats of IEEE 
standard 7S4„ not© that Jbo is explicitly 
included in the string* Interpreting 
the string as a binary number with the 
binary point placed between Jbo and bi, 
th© numerical value of the significand 
is 5, Where 0 < 5 :S 2 - 2^"^. 
Likewise, for convenience it is possible to speak of 
a triplet of miasks corresponding to the above triplet of 
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nimbers. For example, it is possible -co prepare ana 
inanipulate a particular mask as a Bl^nifxcaxid mssk. 

Note that the above splitting of a floating-point 
nxxmbex (or a mask) is for convenience and is 
non-limiting. In particular, it is still possible to 
represent a complete floating-point machine number as a 
single seQuence o± bits, and it is still possible to 
speak of a single mask corresponding to sucH a complete 
floating-point machine number. 

The value which corresponds to such a triplet of 
bit strings is given by: 

1. If e = a^^-l and S ^ 1^ then v is NaN regardless 
of S- 

2. If e = 2^-1 and S = 1, then v « (-1)^ x oo 
5 (Infinity) . 

3. If 0 <: a < 2^ - 1 and S > 1, rhen v = (- 
IJ^ X 2* - X s (Normalised numbers) . 

4. If e = 0 and s < ± tuen v « (-i)** x 2^^^ x s 
(Denormalized numbers and zeroes) . 

Machine numbers are herein denoted in underlined 
italics (such as a machine number a) . 
Rounding 

Mathematically, most numbers in the set of real 
numbers cannot be represented by a finite number of 
digits (not even the entire subset of rational numbers), 
and most of those rational numbers which can be 
represented by a finite number of digits cannot be 
represented by the small fixed number of digits in the 
various floating-point formats. It is common, therefore, 
that the results of a floating-point operation be 
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adjusred^ or "'rounded" to fit within the confines of the 
floating-point representation. The result o£ a rounding 
operation is a floating-point number that approximates 
the precis© value thar should resulr from the 
computation. Rounding, when applied/ always introduces an 
error into tha compuration, bur this error is small 
enough to be ignored in the vast majority of practical 
applications- 

In practice/ the precise value is not necessarily 
computed/ but rather an intermediate result that 
represents floating-point numbers using a larger nizmber 
of significand bits than the permitted format o£ the 
output (also referred to herein as the ^outpur result") . 
The rounding thus consists o£ a truncation o£ the excess 
digits of the intermediate result, and a possible 
incrementing of the least significant bit of the 
significand/ with a possible carry of this incrementing 
to more significant bits of the significand. 

There are several different rounding modes, all of 
which must be taken into account when generating test 
cases r 

• Tound up — round toward the closest 
floaring-point number to the inrermediate result 
that lies between the intermediate result and +00. 

• round down — round toward the closest 
floating-point number to the intermediate result 
that lies between the intermediate result and -00. 

• round to zero — round toward rhe closest 
floating-point number to the intermediate result 
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nhar lies berween the intermediate result and 
zero . 

• round to nearest/even - round toward the floating 
point number closest to the intermediare result 
regardless of the direction. If the intermediare 
result: lies axacrly halfway between the nearest 
larger floating-point number and the nearest 
smaller floating-point: number^ choose the nearest 
floating-point number whose least significant 
significand bit is zero T'even") - 

Coverage 

How does one know that a certain set of tests is 
sufficient? This question is related ro t2ie norion of 
coverage/ that is, to the comprehensiveness of the set 
related to the target floating-point unir- Coverage 
models are usually defined^ and the set of tests should 
fulfill all the existing requirements. A coverage model 
constitutes a set of related cases- Coverage modeling is 
discussed in ^'Software negligence and testing coverage" 
by C. Kaner, Proceedings^ of STAJR the F±£th 

Intex'nat±ansiX con£&retnc&^ Sofzware T&stlngj. Analysis and 
R&vlow, pages 299-327/ June 1996; and ^"User defined 
coverage - a tool supporred methodology for design 
verification'' by R.Grinwald, E. Harel, M. Orgad, S. Ur, 
and A. 2iv, ProcGedings of the 35th Design Automation 
Conference (DAC) , pages 158^153, June 1998. 

As an example, a common coverage model ^ albeit one 
that is far from trivial to fulfill — requir-es 
enumerating all major IBEE Floating Point types 
simultaneously for all opor^Lnds of all FP instrucrions . 
For a given instruction with three operands, say ADD, 
IL 9^2001^-0023 12 



rhis pcitenrially yields in the order of a thousand (10^) 
cases ro cover, assuming 10 major PP types (+NaN'S/ 
±Infinity, iSero^ ±Denorinal/ ±Normal) . This modsl can be 
further refined by adding more FP types / such as Minimum 
and Maximxim Denormals, and so forth- Obviously, nor all 
cases are possible (for example, the addition of two 
positive denormal numbers cannot reach infinity) , so rhat 
the actual number of cases is in fact lower than the size 
of the Cartesian product, A coverage models or rhe set of 
all coverage models, is really an attempt to partition 
the set of all calculation cases in such away that the 
probability distribution should be similar for all 
subsets, 

A Generalized Test-Case Generator 

Consider an auromaric test: generaror whose inpur is 
the description of a coverage models and whose output is 
a set of tesrs covering that model. K coverage model is 
defined by specifying a set of different constraints to 
be fulfilled, where each constraint corresponds to a 
particular task targeted by the coverage model- More 
precisely, a coverage model will have the form of a 
sequence of FP instructions, with sets of constraints on 
rhe input operand (s), the intermediate result (s), and the 
results of the participating instructions. Covering rhe 
mroxtel their ^r^Qn3±re-s — p-rovidtng — a — — of- -tesirs — wh ich 
display the instruction sequence, and which possesses the 
property that each constraint is satisfied by at least 
one test of the set. The general appearance of a single 
instruction constraint is of the following form: 



13 



FPinst (Opl in Pattiarnl) (Op2 in Fatt ern2) (IntRes in 
Patterns) (Res in Pattern4) (1) 



where FPinsr is a generic ^loaring point instruction 
with two inpi3t. operands (Opl and Op2), one inrermediate 
result (in'CRes) / and one output result (Res) . The case of 
two input operands and a single intermediate result is 
used here for simplicity, but of course generalization to 
any number of such parameters is possible. 

A Pattern is a construct representing the logical 
union (u) among sets of TP numbers. The sets serve as 
constraints defining (in fact limiting) the allowable PP 
nxambers for each term of expression (1) . Patterns have 
the general following form: 

Pattern = Seti u Set2 . . . u Set^r (2) 

where each Seti is a set of FP numbers. Each task of 
the coverage model corresponds to a specific selection of 
Seti for each Pattern. Covering the task reduces then ro 
selecr a data -tuple where each individual datum belongs 
to the corresponding selected Set^. Thus, the niomber of 
different tasks originated from a single such instruction 
is the product of the number of Sets for each 
participating Pattern, The number of tasks for a sequence 
~xs the product o£ the number of tasks for each~a:na±vxaual"~ 
instruction - 

There are different kinds of constraints on FP 
ntombers : 

• Ranges — constraints on the upper and lower 
limits of a machine number or element of number 
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rriplet (for example ^ an exponent between 0 and 
2) ; 

• Weights — constraints on the limits of the number 
of bits of value 1 wirhin a machine nximbcr or 
element of a number triplet (for example: at 
least 1 bit ser in the signif icand) j 

• Run-length — constraints on the lengths of 
conrinuous streams of I's and/or O's (for 
example: a stream of at least 45 consecutive I's 
in the signif icand) ; and 

• Masks ~ constraints on individual bits of a 
machine nvimber (described herein in detail) , 

It is also possible to specify a set for which the 
selected value should be a function of the value 
previously selected for another input operand. For 
example, a selected exponent can be at a distance of at 
most 2 from the exponent selected for the previous input 
operand* Set operations (intersection, union, complement, 
of same and different set types) are also possible. 

For a generalised test-case generator ^ any 
architecture resource which might influence FP 
instruction' s results is settable as an input. For 
example, in the non-limiting case of IEEE standard 
architecture/ this applies to Rounding Modes and Enabled 
JE!ljagsL* 

A generalized test-case generator solves constraints 
that are derived from set restrictions on instruction 
operands. Given a restriction^ the generator ideally 
seeks a random instance that solves the restriction, 
where the solutions are uniformly distributed among the 
set of all solutions. In practice, the complexity 
XL 9-2001-0023 15 



involved is sometimes overwhelming^ especially for 
complex or multiple restrictions. In such cases r th© 
generator at least ensures that each solution has a 
reasonable probability of being selected. As described 
above/' constraints can be given on the input operands/ 
the output or even on both simulraneously- ir should be 
clear that there is a significant leap in complaxity 
involved in solving consrrainrs on outpurs- Indeed^ in 
contrast to the case of the input operands^ the 
constraint on outputs includes the instruction semantics. 
However, even output constraints are usually solvable 
analyrically in reasonable time complexity. Constraint 
restrictions start to become largely intractable wh^n 
simulraneous constraints are requested on both input 
operands and outputs. For example, it is largely unclear 
how zo rind an instance in which the result of a MOb 
instruction has at least 45 bits set in its significand 
and the inputs are restricted by specific ranges or 
masks. Such a case might seem artificial, but it is often 
the case that cases such as this one are important to 
check due to specific implementation methods. Moreover, 
during the implementation itself, it is sometimes 
important ro escplor© whether some cases are possible at 
all — it is desirable to know if a solution does not 

exist^^^ Knowing that some cases can be n eglected can 

critical in optimising the microarchitecture. In fact, in 
many cases, ir can be shown that the constraint problem 
is NP-Hard. Thus, the generalised test case generator's 
approach for these problems should be heuristic, mixing 
probabilistic, search space, and semi-analytic 
algorithms. Some Important cases of simultaneous 
IL 9-2001-0023 16 



constraints, howerer, are solvable analytically/ 
including, for eKample, Range constraints or Mask 
constraints on all operands for the FP ADD instruction. 
Test Generation via Masks 

The present: inventors have realized that it is 
possible to advantageously specify rbe generarion of test 
cases using masks to constrain the floating-point numbers 
of the test cases, as described previously. Masks are an 
important way of defining sets of floating-point machine 
numbers by providing constraints on the bits of those 
numbers • 

There are a number of important advantages in 
utilizing masks for specifying constraints in the 
generating of resr cases: 

• Masks are easy to visualise^ understand^ and 
manipulate, using masks does not require the user 
to develop any complicated algorithms or to write 
any software code, as is required for prior-art 
generation functions. 

• Masks can be easily created to describe a wide 
range of machine numbers having specific 
properties of interest to floating-point unit 
designers. It is common for an implementation to 
treat a specific bit, or set of specific bits, in 
a particular manner, and masks are a natural 
means to introduce a bias towards numbers having 
specific bits set to prescribed values, while the 
values of other bits are assigned randomly. 

• At the architectural level, masks allow defining 

all the IEEE generic types of FP numbers and 
symbolic objects: Normals, Denormals, infinities, 
IL 9-2001-0023 17 



zeroes, NaN's, and so forth- Thus, ex coverage 
model handling all types of moinbers and symbols 
is expressible through masks • 
• Although masks ax© not as general as rhe range / 
weight, or run-length constraints previously 
described, ir is relatively straightforward to 
produce a finite {though possibly large) inask set: 
which incorporates and expresses all the 
properties o£ a given range ^ weighty or 
run-length constraint. 
Figure 1 illustrates the development of 
mask-constrained rest cases. (Note that iri this 
particular illustration, a f loating-'point number is taken 
as being associated with a single masik/ as opposed to a 
triplet of masks as previously discussed.) Starting with 
a test concept 101, a specified floating-point operation 
along with some Testing Knowledge 102 is compiled into an 
operation^ mask, and rounding specification 103, which 
selecrs a rounding mode for the given floaring-point 
operation and a set of suitable masks 104, which 
includes r 

• Mask(s) for input operand (s) ; and 

• Mask for the output result. 

Floating-point operation, rounding mode, and mask 
set 104 are input to a floating-point test case generator 
105. Floating-point test case generator 105 in turn 
outputs the floating-point operation, rounding mode, and 
a set of machine numbers lOS, which includes: 

• Machine number (s) for the input operand (s); and 

• Machine number for the output result . 
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Floating-point operarionx rounding mode, and machine 
nuinber set 106 constitute a solution 107. If^ however, 
there exists no set of machine numbers 106 compatible 
with mask set loa {which is possible, such as when the 
specified result mask is incompatible with the input 
operand masks given the specified floating-point 
operation) , then there is no solution. In this case, 
solution 107 is a determinarion that no solution actually 
exists. 

If there exists a set 106, the machine niMbers 
thereof, along with the floating-point operation and 
rounding mode, would then be input to the target 
floating-point unit (not shown) to see if how the unit 
performs for the given floating-point .operation on the 
given inputs, and if the unit produces the given output 
result. If rhe rarger floating-point unit properly 
duplicates the maichine number for the output result as 
given in set 106, then the unit has passed this 
particular test. Otherwise, if the target unit does not 
properly duplicate the machine number for the output 
result as given in set: 106, there is a design error in 
the unit which must be corrected. Ijikewise, if solution 
107 determines that rhere is in fact no solution, it may 
be possible to test the target unit to verify that the 
unit in f.act does not perform any computation that 
results in a set of machine numbers corresponding to the 
masks given in ser 106, because doing so also indicates a 
design error in the unit. Thorough testing requires 
exercising the target floating-point unir with a large 
nxamber of suoh test cases for a variety of conditions - 
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As an sicample (toasiaci on me conceprs in IEEE 
standard 754), consider a hypothetical binary floating 
point format of eight bits, whose structure is seeeffff^ 
Namely, there is one bit for a sign, three bits for a 
biased exponent and four bits for a fraction. In analogy 
with the IEEE formats single and douhl&r "the significand 
has five bits, EjnXx, = -2, E^l^ = bias - 3. Given three 
masks = 0100^101, M^, = OOlj^rlOll, and JVT^ = OlOxxlOx, the 
solution requires three floating point numbers a, h, and 
a, which are compatible with the respective masks, such 
that c ^ round (a + b) . Assuming that round stands for 
round to nearest /even, one solution is a = 01000101, 
b = 00101011, and c = 01001100, 

In the sch^e illuatrared in Figure 1, it: is 
necessary to construct floating-point test-case generator 
105 with the following properties in mind: 

• All valid solutions must have roughly the same 
probability of being produced by the test-case 
generator; and in particular, 

• The solution set must not arbitrarily exclude any 
valid solution. In other words, every valid 
solution must have a non-zero probability of 
being selected. 

Currently, however, there are no such 
mask-constrained floating-point test-case generators 
available, even for a restricted set of floating-point 
operations such as addirion and subtraction. In order to 
create a f loating-'point test-case generator, it is 
necessary ro have as a minimum, a floating-point 
test-case generator for addition and subtraction, which 
can solve the following problem: Given masks for three 
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machine numbers and a rounding mode, generate machine 
numbers a, h, and which are compatible with the given 
masks and satisfy a = roundly ± ^) , where round 
corresponds to the given rounding mode. 

There is rhus a need for, and it would be highly 
advantageous to have, a mask-constrained f loaring-poinr 
test case generator for floating-point addition and 
floating-point subtraction which has rhe desired 
properties listed above- This goal is met by the. present 
invention • 

SXJMMRRY OF TBS XKVENTXC3ISr 

The present inventors have recognised that the lack 
of a practical framework for generating constrained, 
meaningful test cases is a major deficiency in the prior 
art, and represents rhe principal obstacle ro efficient 
verification. The present invention, therefore, 
approaches this problem from the overview presented in 
the following sections. 

It is an object of the present Invention to solve 
the following problem: Given masks for three machine 
numbers and a rounding mode round, generate machine 
numbers, a, b, c, which are compatible with the maslcs and 
satisfy a = roui^dia ± h) , 

It is an object of the present invention to develop 
a method for generating three floating-point machine 
numbers a, and a corresponding to three given masks 3^, 
Mb, and Mo, and a rounding mode round, such that 
c = round (a ^ b) , where either or both a and b may be 
posirive or non-positive - 
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Ir ±s also an object of the present invention that 
all valid solutions have roughly th© same probability of 
being produced by the method^ and that no valid solutions 
be excluded • 

It is moreover an object of the present invention 
that the method support: general binary f loarlng-point 
standards, including^ but not limited to IEEE standard 
f loating-poinr aritiimetic when generalized to include all 
allowed PP format sizes (such ae 32 bits^ 64 bits, 80 
bits, 128 bits) . With respect to this particular goal/ it 
is noted that the innovations of the present invention 
are not necessarily limited to binary implementations/ 
but can be applied in other niixnber systoms as well. 
Although binary systems are of primary importance/ 
because current floating-point standards are expressed 
explicitly in terms of binary number©/ it is understood 
that the present invention is not limited to binary 
arithmetic/ and the details of the embodiments herein 
presented are to be taken as non-limiting examples. 

It is further an object of the present invention to 
develop a system for implementing tne method- 

The basic problem may be split into two 
sub-problems, which may be solved by rwo generators of 
machine nuinbers/ respectively: 

• Ploating-point generator for addition : Given six 
masks, for biased exponents and for signif icandS/ 
Msa, Met>, Mst,, Meat Mscg the generator either 
generates three non-negative machin<y numbers, a/ 
b/ C/ whoso biased exponents and whose 
significands are compatible with the 
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corresponding masks and satisfy/ 

c - round + , or states that there is no 
i solution. 

• Float inq~polnr generator for subtraction ; Given 
siK maskS/ for biased exponents and for 
significands, M^^, Ms^, M^h, Msh, Mecr Msa the 
generator either generates rJixee non-negative 
machine numbers, «/ b, c, whose biased exponents 
and whose significands are compatible with the 
corresponding masks and satisfy, 

c = rovnd(a - b) , or states that there is no 
solution. 

Not© that when the problem is reduced to the two 
cases above^ the three machine numbers, a, b, and c, are 
all non-negative. In the case of the floating-point 
generator for subtraction, the values assigned to a and b 
are interchanged if necessary so that a ^ Note also 
that in the case of the floating-point generator for 
subtraction, where a 96 ^ it is possible to exercise the 
target floating-point unit with two test cases resulting 
from a single solution found by the generators one test 
case is simply c = round (a - b) , where c is positive; and 
the other test case is c = round ib - a), where c is 
negative • As far as the floating-point test-case 
generator is concerned, these are mathematically 
identical and do not require a separate method Of • 
solution. From the standpoint of verification testing of 
the target floating-point unit, however, they are 
distinct problems, and the successful passing of one test 
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by rhe target f loa-ting-point unit does not necessarily 
imply the successful passing of the other test. 

The method of the present invention solves the above 
problem fox simulraneous mask constraints. That there 
is a mask constraint on both input operands £ and Jb, and 
on ttie ourpur result c for any specified rounding mode 
and for the floating-point operations FAPD and FSUB. It 
is noted that the distinction between addition and 
subtraction can be implied in the sign of the operand bz 
If b < 0, a STobtraction operation is implied; otherwise , 
an addition operation is implied. It is further noted 
that a single solution found by the binary floating-point 
test-case generator of the present invention can be 
verified on the target floating-point unit for both SADD 
and S-SOB by adjusting the sign of one of the input 
operands- That is, given the solution a, b, and a where 
a = roundffa ^ b) , it is possible to verify FADD on the 
floating-point unit for a, b, and and it is also 

possible to verify fsot on the floating-point unit for a, 
-b, and £• The internal implementation of FADD in the 
target floating-point unit can be structurally different 
from that of FSUB, justifying a separate verification 
pass, but the identical result c should be output in both 
cases . 

Rounding Modes Meeded 

In the general case, there are the following 
rounding modes; 

roujad € (round down, round up, round toward zero, round 
to jaeare^t/evan} (3) 
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However^ since rhe problem has been reduced to cases 
where the three machine numbers are all non-negative # it 
is possible to omit, without loss of generality, the 
round toward z^ra mode- Since all non-zero machine 
numbers will be positive, this mode is equivalent to 
round dotnrn. 
Novelty 

There are several novel aspects of the present 
invention : 

• The present invention goes beyond the employment 
of masks for expressing floaring-poinr machine 
number constraints to include a novel use o£ 
masks for controlling rhe stream of carry bits 
created' during the floating-point addition 
operation- In addition to providing a key portion 
c£ the method/ this also enriches the 
verification prac&ss^ by allowing a floating-point 
unit to be exercised through different carry 
configurations . 

• The present invention discloses a fixed-point 
generator which has a number of unexpected uses 
in various places throughout the floating-point 
rest-case generation. 

Notational Conventions 

The following notation is employed herein in the 
context c = round [a ± b) : 

a, Jb^ a floating-point machine numbers 

ei biased Gxrpon^n-t 

Si unbiased escponent 

M± mask 
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n boundary indeX/ the position of the leftmost 2 

in a carry sequence 
N number of bits in a fixed'-point machine nuinber 
jp number of bits In a slgnificand (including the 

most significant bit) 
q maximum definire binary point shift, usually 

comparable in magnitude to p (shifts greater 

than q are considered Indefinite) 
q± biased exponent shift — 

e^y Qi> = So - (1) 

Qi unbiased exponent shift - 

Dd. = JSTc - jSs, Qb - Ea - Ejt, (2) 

S± signif icand 

w number of bits in a biased exponent 
y, s fixed-point machine numbers 

Figure 2 is a block diagram of a floating-point 
addition/subtraction test-oase generator according to the 
present invention. A mask set 201 is input into a 
processing unit 202, which generates test-cases in an 
output 250^ which, as previously noted, may be either 
floating-point machine numbers compatible with mask set 
201 and a specified rounding mode (not shown) or the case 
of ''^no solution". In the non-limiring example presented 
herein, processing unit 202 includes a g-g, g^ / 
selector 205 whose output goes ro a biased exponent 
generator 210, which includes a definite biased exponent 
generator 215 and an indefinite biased exponent generator 
220. For addition, g^, / q^, q^ selector 205 selector 
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selecrs and q^, whereas Tor subtraction tgr^, 05? / qc 
selector 205 selector selects q^, and q^ Then, a Q^, Qt, / 
Qh, Qa calculator 225 computes Qa and C2i> for addition or Qi^ 
and for subtraction, and sends rhem to a signlficand 
generator 230, which includes an addition significand 
generator 235 and a subtraction significand generator 
240. Both biased exponent generator 210 and significand 
generator 230 rely on a fixed-point generator 245 in 
making their respective computations. Finally, the output 
of biased exponent generator 210 and significand 
generator 230 are combined to output solution set 250, 

Nore that/ in Figure 2, biased exponent generator 
210 is a non-limiting example of a more general exponent 
generator thar need not be restricted to biased exponents 
only. Because an unbiased exponent shift is calculated, 
it Is also possible to conceptualize the system o£ 
Figure 2 in terms of such general exponent generator. 

Given the selected and calculated values for q^f qt^ 
Qa^ Qb for addition (or q^f qc, Qhf £?c/ in the case of 
subtraction) , it is possible to produce the biased 
exponents and rhe significands independent ly, by invoking 
number generators such as those shown in Figure 2, and 
defined as follows: 

Definite biased exponent generator 215 performs the 
following: Givan the two non-negative integers gi, q^^ 
with qi 6 {0, 1}, and three masks of length w. Mi, Mzi M3, 
for the biased exponents, definite biased exponent 
generator 215 either generates and outputs three biased 
exponents ©1, e^r ©3/ which are compatible with the 
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respective masks and sat;isfy es ei + <3*i = ©2 + <3r2f or 
outputs that no solution exists • 

Indefinite biased exponent generator 220 performs 

the following: Giv^n thia two non-n^^gativa integers q^^ 

With gi € {0, 1}, and three masks of length Wf Mif Ms, 
for the biased exponents/ indefinite biased exponent 
generator 220 eirher generates an integer q[2 and three 
biased exponents ei, e2/ g^, which are compatible with the 
respecrtive masks and sarisfy ©3 = 61^-31 = 65 + Q20 such 
that qz > q, or outputs that no solution exists- 

Addition significand generator 235 performs the 
following: Given two non-negative integers Qa, Qt, one of 
which is 0 or 1/ three masks of length p for the 
signif icands, Msat Mst>f ^scr a rounding modo^ 

roimdf € (rouiid dawn^ jcoundt up^ round zo nearest/evezi} f 
addition significand generator 235 either generates three 
signifioands Sor Si,, 5c which ar© coznpatil^le wirh the 
respective masks and satisfy 5c = round (2'^^ 3^ + 2"^ 3^) , 
or outputs that no solution exists* 

Subtraction significand generator 240 performs the 
following; Given two non-negative integers a, Qc, one of 
which is 0 or 1, and three masks of length p for the 
signif icands, Ms^f Mshf Msesf and a rounding mode, 
round e {round dovirnj, round up, round to nearest/even}, 
subtraction significand generator 240 either generates 
three signifioands $3, $h, Sc, Which are compatible with 
the respective masks and satisfy 2'°°'Se = round (Sa - 2" 
Sh) , or outputs that no solution exists. 

Figure 3 is a simplified block diagram of 
fixed-point generator 245 (as also appears in Figure 2) 
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according -to the present invention. A mask set 303 
containing laasks M^, Myf and Mg^ for three f isced-point 
machine n-Qinbers y, and z is input into fixed-point 

generator 245. In addition, a mask 305 iTor rhe carry 
Stream is also input. Note that M^, My, and Hz each have N 
characters corresponding to th© i\7 bits of and z, 

whereas Mc has N + 1 characters* Based on the inputs, 
fixed-point generator 245 outputs a solution 307^ which 
either contains a set of three fixed-point machine 
numbers y, and si which are compatible with mask set 
303 {Mjc, My, and jS4) and mask 305 (Mc) ; or which is the 
case of ^^no solution-^. 

It is to be noted that all of the generators of the 
present inxrention generare a solution rhat is in effect 
randomly selected from the complete set of valid 
solutions, and if rhe ser of valid solutions is empty, 
the case of ^^no solution'' is output instead. It is not 
required that a generator output all of the valid 
solutions, but only a single solution, provided that at 
least one valid solution exists, it is also not required 
that a generator be able to determine how many valid 
solutions there are. The generators / however, do not 
arbitrarily exclude any valid solution. That is, there is 
a non-zero probability that any specific valid solution 
will be generated. Subsequent operation of a generator 
will usually generate a different, randomly-generated 
solution. If the complete ser of valid solutions is of 
high order, then it will be overwhelmingly probable that 
repeated operations of a generaror will result in the 
output of distinct solutions, so repetitions of identical 
solutions are normally rare. 
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Solution-SeeTcing Merhod 

As will be discussed in de-bail below^ the generators 
make extensive use of a solution-seeking method 
illustrated in Figure 4- In the present invention^ the 
5 solution-seeking znethod i$ used in a nuinber of places to 
seek a solution that may satisfy any ot a variety of 
specified mathematical conditions. The solution that is 
O sought may represent intermediate results wirhin 

m processing unit 202 (Figure 2) which are used in 

M 10 constructing floating-point test-case solution 250^ as 

^ well as floating-point test-case solution 250 itself. 

As iliusrra-ceci in Figure 4, associated with the 

■i I 

!s1 method is a finite list of input choices 401 which 

ry contains a set of parameters for the solution being 

'j; 15 sought. Other parameters may be applied, as discussed 

flj below. The method begins at a start block 403, 

immediately after which a decision point 405 checks to 
see if list 401 is empty. If so^ then there are no 
available parameters for the solution, indicating that 
20 there is no solution/ in which case at an output block 
407 the method terminates with the output of the ^'no 
solution'' case. Otherwise^ in a selection block 409 one 
of rhe elements of list 401 is selected at random. Then^ 
in a solution-searching block 411 a solution is sought 
25 utilizing the selection from the list. A decision point 
413 checks to see if a valid solution has been found. If 
so, the solution is output in a block 415^ and the method 
terminates. Otherwise^ in a block 417 the selection made 
in block 409 is erased from list 401^ and the method 
30 resumes once again at decision point 405. Ultimately, 
because list 401 is finite, either a valid solution will 
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be re-turn^d in block 4lS, or the '*no solution" case will 
be returned in block 407 • 

The following points are noted regarding 
solution-seeking block 4X1: 

• The inner details of block 411 depend on the 
specific derails of the problem for which a 
solution is sought; 

• A valid solution may depend on additional 
parameters within block 411 besides the selection 
made from list 401; and 

• Block 411 may encapsulate one or more instances 
of this same solution-seeking method — that is, 
the entire solution-seeking method may be used 
recursively within block 411, such that another 
set of all the elements shown in Figure 4 appear 
within block 411. For example, the floating-point 
addition test-case generator as shown in 
Figure 2 utilizes a solution-seeker as shown in 
Figure 4, whose list 401 contains a finite number 
of qa, qb pairs. But within the floating-point 
addition test-ease generator rhere is addirion 
significand generator 235 which itself utilises a 
solution-seeker as shown in Figure 4, whose list 
401 contains a finite number of elements (^^tails 
triplets") which are based on a Qa/ Qh pair 
calculated from the q^, pair selected from the 
f loaring-poinr addirion test-case generator's 
list 401 • This is explained in detail below^ but 
is mentioned here to show how the solution-seeker 
of Figure 4 may be used recursively. 
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• Block 411 and decision-point 413 may be repeared 
seQuentially within the loop shown in Figure 4f 
with found solutions from one block sent to the 
nent such block as an intermediate or partial 
solution. 
£;K;aTnple of a Solution 

Herd XB an example of a triplet of masks which may 
be used as input to a binary floating-point addition 
test-case generator according to the present invention, 
illustrated in the non-limiting case of IEEE standard 754 
double precision (1 sign bit/ 11 biased exponent birs/ 52 
fraction bits, for a total of 64 bits) ; 

The furnished masks specify two denormal floating 
point niimbers as input operands for addition, and the 
mask for the sum specifies a normalized floating-point 
ntimber. 





3 




3 




0 


00000000 
000 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxx 


Mb 


0 


00000000 

OOP 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxx 




0 


00000000 
001 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxx 
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One possible solu-blon Is: 





s 


e 


s 




0 


00000000 
000 


1111111111111111111111111111111111111111 
111111111111 


b 


0 


00000000 
000 


0000000000000000000000000000000000000000 
000000000001 


a 


0 


00000000 
001 


000000 00 OOOOOOOOOOOOOOODOOOOOOOOOOOOOOOO 
000000000000 



Another solution is: 





s 


e 


s 


a 


0 


00000000 
000 


1111100000000000000000000000000000000000 
000000000000 


b 


0 


00000000 
000 


0001010101010101010101010101010101010101 
010101010101 


c 

. 


0 


00000000 
001 


0000110101010101010101010101010101010101 
010101010101 



The complexity of zhe method according to the 
present invention is polynomial, and tests confirm that 
it operares efflcianrly and quickly: in practice^ 
solutions are found almost immediately for a large 
v^rUiety Of' ins1:a"nGesT ~ " ~ ~- - ~- 

It will also be understood that the system according 
to the invention may be a suitably programmed compurer- 
Thus, the invention contemplates a computer program being 
readable by a computer for e2s:ecuting the merhod of the 
invention. The invention further contemplates a 
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machine-readable memory tangibly embodying a program of 
instructions in data storage executable by a machine for 
performing the method of the invention. The invention 
moreover contemplates a machine-readlable memory tangibly 
embodying a program o£ instructions in data storage 
executable by a machine for emulating the system of the 
invention. 

According to the present invention there ±s provided 
a system for generating floating-point test-cases for 
verifying the operation of a floating-point arithmetic 
unit, the system including a processing unit which 
Includes: (a) an exponent generator^ for generating 
floating-point exponents; (b) a significand generator, 
for generating floating-point isignif icands; and (c) a 
fixed-point generator coupled to the eacponent generator 
and to the signficand generator,- wherein the processing 
unit is configured to receive a specified arithmetic 
operation selected from a group that includes addition 
and subtraction, a specified rounding mode, a first input 
operand mask, a second input operand mask, and an output 
result mask; and wherein the processing unit is 
configured to output a set of floating-point numbers 
which includes a first input operand compatible with the 
first input operand mask, a second input operand 
GompatJ-ble_ with _the second^ _input _ operand, mask/^ _ and_ an 
output result compatible with the output result mask; and 
wherein the output result corresponds to the specified 
arithmetic operation on the first input operand and the 
second input operand for the specified rounding mode'. 

In addition^ according to the present: invention 
there is provided a method of seeking a solution, if a 
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solution exists, to a specified mathematical condition, 
whertjin the solution is used in constructing a 
floating-point test-case for verifying the operation of a 
floating-point arithmetic unit, wherein a complete 
generated test case ±s a set of floating-point nioiubers 
for a specified arithmetic operation selected from a 
group including addition and subtraction, and for a 
specified rounding mode, and wherein a generated test 
caase includes a first input operand, a second input 
operand, and an output result; and wherein the first 
input operand is compatible with a first input operand 
mask/ the second Input operand is compatible with a 
second input operand mask^ and the output result is 
compatible wirh an output result mask; the method 
including the steps of: (a) preparing a list of choices 
upon which the solution is based; Cb) testing whether the 
list of choices is empty; (c) outputting^ i£ the list of 
choices is empty ^ that no solution exists; (d) randomly 
choosing/ if the list of choices is not empty, a choice 
of the list as a selection; (e) searching for a solution 
to the specified mathematical condition, based on the 
selection; (f) outputting, if the searching was 
successful, the solution; (g) erasing, if the searching 
was not successful, the selection from the list; and (h) 
_ repeating, step_ (^)_ tJajQugi _:Stap_ (g)__ .until joutput±lng_ 
occurs . 

Furthe2rmor&, according to rhe present invention 
there is provided a method of generating a set of 
fixed-point numbers containing a firsr addend, a second 
addend, and a sum, wherein the first addend is compatible 
with a first addend mask, the second addend is compatible 
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wirh a second addend masky the siom is compatibl© with a 
sum mask, and wherein the addition of the first addend 
and the second addend results in a carry sequence of 
carry bits, wherein eacn carry bit has a unique index in 
5 the carry sequence, wherein the^ carry sequence is 
comparible with a carry sequence mask and wherein each 
ri carry bit has a value in th^ group consisting of 0^ 1, 
O and 2r and wherein rhere exists a boundary index in the 
fn carry sequence corresponding to the lowest index of a 
1^ 10 carry bit having the value 2; the method including the 
n steps of: (a) constructing a lisr of possible boundary 
indices; (b) testing whether the list is empty ^ (o) 
^ outputtlng, if the lisr is empty, . thar no solution 
fU exists; (d) randomly choosing, if the list is not empty^ 
!^ 15 a boundary index from rhe list as a selection; (e) 
fy searching for a carry sequence based on the selection, 
which is ccmpatible with the carry sequence mask; (f) 
erasing, if the searching was not successful, the 
selecrion from the lisr; (g) const rucring, if the 
20 searching was successful, a first addend compatible with 
the first addend mask, a second addend compatible with 
the second addend mask, and a sum compatible with the sum 
mask; (h) ourpurring the first addend^ the second addend, 
the sum, and the carry sequence; and (i) repeating step 
~ - - -(sl) .through -St ep^ (hj until^autputting ocaurs.- _ ^ _ _ 

BRISF DBSCKZFTZON O^" OIHB BR2^S7X1T6S 

The invention is herein described, by way of example 
30 only, with reference to rhe accompanying drawings ^ 
wherein: 
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Figure 1 illustrates the progression from test 
concept to mask specification to the finding of a 
solution^ or the determination that there is no solution. 

Figure 2 is a block diagram of a floating-point 
addlrion/subt:raction test-case generator according to the 
present invention . 

Figure 3 is a simplified block diagram of a 
fixed-point generator according to the present invention- 

Figure 4 is a flow-chart showing the steps of a 
generalized solution-seeking method according to the 
present invention. 

Figure 5 graphically shows rhe finding of a solution 
for the indefinite biased exponent generator* 

DB0CRXP3>XQN OF TTOB PX(SRBXVEffiD SMBOPXMSN*FS 

The principles and operation o£ a floating-point 
addition/subtraction test-case generating method and 
generator according to the present invention roay be 
understood with reference to the drawings and the 
accompanying description. 
Outline of rhe Method 

As was stated earlier^ the problem of generating 
- floating -pointr- nianbers; ~sati-sfying^ "c"- roui^cf(&~± b) ,' may" 
be divided into two cases: addition of non-negative 
machine numbers and subtraction of non-negative machine 
numbers . 

Consider first the addition case. Namely, let 
c « roundlB + Jb) , where a, Jb, c are non-negative machine 
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numbers* The biased exponent shifts are denoted ^STa = ~ 
©a snd = ©c - &£>^ where e^r are the biased 

exponents. It is not difficult to see that g^, gi> are 
non-negarive integers, one of which must: be either 
0 or 1, Likewise/ the unbiased exponent shifts are 
denoted ^ Se ~ E:^ and Qja = £o - Eh, where Ea/ JSi, and Ec 
are the unbiased exponents- It is easy to see that Oa, Qb 
are also non-negative integers, one of which must be 
either 0 or !• Usually = q„, gj, = bur this is not 
always so. For addition there are the following five 
cases: 

> 0, > 0, ee > 0: Qa ^ Qa/ Qh =^ Qh (S) 
©a = 0,&i, > Q, &a > 0: Qa = qfc - 1/ Qh ^ qb H) 

> 0, ei, = 0,eo > 0: ^ q^, « - 1 (8) 
= O^ej, = 0,ec > 0: = g„ - 1, Qi> <Zb - 1 (9) 

©a = O^Sb = 0,ec = 0: = = (10) 

For subtraction the outcome is similar: With 
c = round i a - there are Qi, = - ei,^ g^ = - e^, 

Qjs = - fJis, and Qc - Ea - £?c/ where one of qj^, q^ is 
either 0 or 1^ and one of Qj^r Qc is either 0 or 1. Thus, 
there are an additional five cases; 

> 0, Si > 0, > 0: ^ q^, « {11) 



O/Si? > 0^ ea > 0: Qe ^ ge - 1/ 
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0, - 0/ ©if > 0: = gc - Ob - 1 



(13) 
(14) 



As illustrated in Figure 2 for the non-limiting 
5 example of processing unit 202, after mask set 201 has 
been given^ the generation of machine numbers is started 
by selector 205, which selects values for qa, (3b for 
addition, or qt,, qc for subtracrion. With known values for 
these numbers, it may be possible for biased exponent 

10 generator 210 ro find valid solutions for the biased 
exponents, after which calculator 225 can calculate Qo, Oi? 
for addition or Qb, Qc for subtraction. For addition, 
calculator 225 relies on Case (6), Case (7)^ Case (S), 
Case (9), or Case {10), above, as appropriate- Likewise, 

15 for subtraction, calculator 225 relies on Case (11), case 
(12), case (13), Case (1), or Case (15), above, as 
appropriate- Based on this, signiflcand generator 230 is 
able to complete the process for output 250. 

The addition/subtraction test-case generator 

20 illustrated in Figure 2 makes use of the solution-seeking 
method illustrated in Figure 4. In order to do this, it 
_ _ is first necessary to produce list 401 for g^, qt in the 
case of addirlon, or q^, qc in the case of subtraction. 
Table 1 shows a list of such pairs: 
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1 


2 


3 


<7 




2 


3 


. . q 








0 


0 


0 


0 


1 


1 


. . 1 


1 



Note thaty while Table 1 ±s expressed as gfj, pairs 

a 

(for addition) / the pairs are equally valid as c[o, gb 
pairs (for subtraction) - Solutions for pairs 0, >g; 

1/ >qp >qt 07 and >q, 1 are handled by indefinite biased 
exponent generator 220 (Figure 2) . Solution for all other 
pairs are handled by definite biased exponent generator 
215. The exact value of q in the table is not specified/ 
but if a convenient value close to p is selected^ there 
will be approximeitely 4p pairs in Table 1- Within the 
f loaring-poinr test-caso generator^ as shown in Figure 2^ 
calculator 225 calculates Q^, Qh for addition or Qtr Qc for 
subtraction^ as may be necassary for internal finding of 
significand solutions^ as discussed below . 

The solution-seeker of Figure 4 is then invoiced, 
with solution-searching block 411 set up to seek 
solutions according to the following (making reference to 
Figure 2) s 

1. use either definite biased exponent generator 215 
or indefinite biased exponent generator 220 as 
appropriate (according to the above criteria for 
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•the inpur pair) ro seek a solution triplet for &af 
e/>/ and Be that is coinpatible with the exponent 
masks of mask set 201; and 
2. oomput© Qstg QjD / Qbr Qa -^isi calculator 225 and use 
either addition significand generator 235 or 
ssubtraction significand generator 240 as 
appropriate to seek a solution triplet for Sg, S^^, 
and Sa rhar is compatible with the significand 
masks of mask set 201. 
If either one or both of the above solutions do not 
exists decision-point 413 (Figure 4) branches to erase 
block 417- Otherwise^ if both of the above solutions 
exist/ the method has found a valid solution triplet e^. 
Sat ^br ^z>; and e^. So, having both exponents and 
significands, and therefore decision-point 413 branches 
to ourpur block 415- 
The Fixed-Point Generator 

As noted previously, several of the computations 
required for a floating-point generator according to the 
present invention require a fixed-point generator. 

In implementing significand generator 230 
(Figure 2) , for example, note that/ since 2^, Cjs {ox Qj^f 
Qa) are known (as output from calculator 225), it is 
possible to shift the significands until they are 
properly aligned so that they have identical exponents, 
and then add (or subtract) them precisely the way 
f isced-point numbers are added (or subtracted) . To do 
this, it is necessary to have a fixed-point generator* 

Before specifying the exact function of the 
fixed-point generator, first consider the process of 
adding two positive binary integers, x + y = The 
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addends in this process are x and y/ and the sum is z. 
Th€: process starts by adding the rightmost (least 
significant) bits of x and y. If the sum is less than 2 
then it is equal to the rightmost bit of 2 an(3 rhere is 
no carry- If the sum is not less than 2, there is a carry 
of 1. Next^ thd carry is added along wirh the following 
bits of X and y. Once again, if the sum is lees than 2, 
there is no carry. If the sum is nor less than 2, there 
is a carry of 1 . This is repeated through the final (most 
significant) bits of x and y. Thus^ during the addition 
process / a sequence of carries/ each of which is either 0 
or 1, is generated- The carry sequence represents the 
carries frcm the successive digits of the addends , Note 
thar rhe fixed-point generator considers only the 
addition of two nuzribers. 

in the following discussion^ the bit-numbering 
convention for a binary fixed-point number containing N 
bits is as follows: the leftmost (most significant) bit 
is assigned an index value of m ^ 0, and the righmosr 
(least significant) bit is assigned an index value of 

If the values of the bits of the addends are Xm — i/ 
^ 3 ayid rhose of rhe sum are Zj/i * "Chen there is the 
equation: 

1 -h J + Cj^i ^ k ^ 2C^ (iZ2 = 0, 1, ...^ i\r-l) (16) 

where Cm is the carry sequence, representing the 
carries resulting from the addition of successive digits 
of the addends . The bits themselves always are such that 
1, jg k € {0^ 1), For the carry, it is also noimally the 
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casts that Cja e {0/ 1}. However, a around up pxaa^ss may 
add an additional 1 to the carry and produce an effective 
carry of 2^ so for this reason it is convenient to allow 
Cffl, Cjo+i s {0, 1^ 2}- Note rhar whereas {m ^ 0, 1, ^ 
1) for the bits of the fixed-point nixmbers, 
(m ^ 0, 1/ N) for the carry sequence. Cm is the carry 

out of bit m and into bit m-lj, whore bit m and bit: jn~l 
actually exist. For example. Cat is the carry into bit ^-1 
only, because there is no bit iV, Likewise, Cq is the carry 
out of bit 0 only^ because there is no bit -1. Co and Cj^ 
are boundary values, and usually both have a value of 
zero. For generality^ however, it is possible that Co ^ 0 
and Cw at 0- 

It is important to note that, while 
Cpi, Cjii+t s {0/ 1/ 2}, there are restrictions on the 
appearance of a value of 2 in a carry position; 

= 2 =^ C^l « 2 (17) 

because the only way a carry position (in m) value 
can be 2 is if the previous carry position (in m+1) is 
Restriction (2) can be rewritten as 

Cj^i < 2 =^ Cm < 2 (18) 
And 



C«4.i = 2 C^ > 1 (19) 

because even if i + j* = 0 in position m, a carry of 
2 in position m+l will propagate at least to a value of 1 
in position in. 
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Because of Restrictions (17)^ (18) ^ and (19)/ it Is 
easy to see that i£ there are any 2's in the carry 
sequence, they are all grouped together to the right, and 
that the^^e is a 1 to the immediate left of th^ loftmost 
2. For example, rhe following is a possible carry 
sequence Into a 16-bit sum: 

0110101222222222 

The following, however, are not possible carry 
sequences; 

0110101222222221- violates 
Restrictions (17) / (18) 

0110100222222222 - violates 
Restriction (19) . 

Thus although Cjn, 0^^^^ e {0, 1, 2} iinplies that there 
are 3 i:- 3^) differenr coiribinations for the pair ic„,C^+i), 
there are in actuality only S different allowable 
combinations / because rhe pairs (0,2), (2,0), and (2,1) 
are ruled out by Restrictions (17), (18), and (19)/ 
leaving only the following set of possible pairs for 

(C„, O^i) {(0,0), (0,1), (1,0), (1,1)^ (1/2), (2,2)} 

(20) 

Set (20) is important in the construction of carry 
sequences. An innovation of the present invention for 
solving the fixed point generator problem is to first 
construct the sequence Cm, and only later construct the 
bits of X, y, and z. This is discussed in detail below. 

The input to the fixed point generator includes 
laasks of length N: m^, My, Mz, of the form described 
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earlier, for the numtoers y, z. The input also includes 
a masky Ma of length iS7 -J- 1, which corresponds to the 
sequence of carries. This carry sequence mask can include 
the characters ^0', ^1'^ ^2', and ^yJ , wnere in analogous 
5 fashion to the previously-defined mask characters, ^0' ^ 
U'/ and ^2' completely specify the value of the 
corresponding carry/ whereas an leaves the 

corresponding carry undetermined. 

The fixed point generator xa then defined as 
10 follows: 

a Fixed-point generator : Given three masks, i^/ My^ 
of length AT for three fixed-point binary 
numbers of ^ bits each/ and given one mask. Mo, of 
lengrh AT + 1 for a corresponding carry sequence, 
15 the fixed point generator either generates rhree 

fixed-point binary numbers k, y, zg which satisfy 
z ^ K y in conjunction with a carry sequence, 
all of which are compatible with their respective 
maskS/ or states that there is no solution. 
20 The operation of the fixed-point generator is 

discussed in the following sections. 

Mask Combination Numbers and Case Numbers 

2S The basic relations which control the construction 

of the sequences = Ym 3§ Zm =^ 

(m ^ 1, , M-l) and c^^ (m = 0 , 1, . . , , N) are the 

condition of compatibility with the masks and Equation 
(16), previously discussed- Clearly rhese conditions 

30 might be self-oontradictory- Where such contradictions 
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exist,, the fixed-point generator states that there is no 
solution. 

Given an index each value of the bits i, j, k 
oorrc5pc>nds to a character in thd appropariate mask* This 
character may be either an ^x' or a number (^0' or ^l')» 
With such a classification of the characters of the mask/ 
each triplet of masks elements is one of eight possible 
types of triplets (for examplei, all three of the 
characters may be ; i corresponds to a number and j 

and k both correspond to an ; and so forth) . Each of 
the eight types of triplets may be assigned a number, 
which is denoted as MCN (Mask Combination Number) , Table 
2 below lists the values. In this table n means a number 
character in the mask fO" or *1') and x means an 
character in the saciask- 
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a5al>Xe 2: Mask Crmlni rta tlcaa IflPaiDibers 

Note that Equation (16) can be alternately expressed 
in the form: 

X j k 2Ca - qa+i (21) 

Given the masks M^r Mr and a numerical value for 

the index m, it is possible to assign an MCN value and 
the numerical values of some of the variables If j, k* 
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For MCN = Dg 1, 2, 3, 4, 5^ and 6^ there is sufficient 
inf oritiattion to compute a Case Number CM, where 



CBT = 1 + J - Jc^ When all rhree 1, J, and ;r are known (n) 

- MCN 0 (22) 

caar = 1 + when only 1 and j are known (n) and k is 

unknown (jc) — MCN J (23) 

CXi - ± ^ k, when only i and k are known (n) and j is 

unknown (x) - MCN 2 (24) 

CKT = 1^ when ooly i is known (a) and Jfe* and j are unknown 

(Jf) - MCN 3 (2S) 

CN = J - -fe"^ When only j and ic are known (n) and i is 

unknown (jc) — MCN 4 (26) 

CN = When only j is known (n) and Jt and 1 are unknown 

(«) - MCN 5 (27) 

CN = -Jf, wnen only is known (n) and 1 and j are unknown 

(x) - MCN 6 (28) 
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Th.n the pairs (C.,C»..) which are possible for 
MCN, oi combination are as follows: 

CW = -1: (0,1) 
^ CN - 0: (0,0-), (1,2) 

CN - 1: (1,1) 

^ ^ 2: (1,0), (2,2) 
• MCM = CW = i + J 

C»= 1: (0,0),(l,l,,(i^2j 
2: (1,0), (1,1), (2,2) 

k 

,5 (0'0)/(0,l),(l,2) 

0= (0'0).(1,1),(1,2) 
^ 1: ^1/0)/ (1,1), (2,2) 

• *«CN « 3, CN = i or MCN = 5, CN = J 
= 0: (0'0)W0,l),(i,i),(i,2) 
"^''^ ^°'°''^^'0)'U,1), (1,2), (2,2) 

• -1: f0'0).(0,l),(l,i),(l,2) 

^-0: (0.0), (1,0), (1,1), (1^2), (2,2) 
• MCN = 7, Ctf^o 

est 

^°'°^'^''°>'^0'lJ'fl.i),(l,2),(2,2) 
arable 3. (C.,C^,) pair3 MOf and CN valties 

The determination of the (c r \ 

j;:;:/;-^^ ^ 
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all three are known, so the various combinations can be 
computed precisely. In the case i = j = 0, k = I, C3a = - 
1. According to Equation (21), the only values of C«, Gn+i 
which can satisfy -1 = 2Ca - Qb+i are Ca - 0 and C^+i = 1, 
5 so the only possible (CC^i) pair is (0,1), as shown in 
Table 3 for C» = -1 under MOS 0. As another example, for 
MCM 1 only 1 and j are known, and k is indeterminate and 
thus CN is simply given by. i + J/ as shown in Equation 
(23) . For i = j = 0, CN = 0, and there are three (Cj,i,Cart-i) 

10 pairs that give this result in Equation (21): (0,0) with 
k -» 0; (0,1) with k = 1; and (1,2) with Jc - 0. These are 
the values shown in Table 3 for cai = 0 under MCN 1. The 
rest of Table 3 is compiled in the same manner. For MCN 
7, a value of CM = 0 is assigned, because all three of i, 

15 j, and k are indeterminate. For this value, the complete 
Set (20) is present. 

This list of 15 CN values in Table 3 is exhaustive, 
because for each pair of MCN, CN all possible pairs (Cm, 
Cart-x) are included for all possible values. This list is a 

20 basis for the construction of feasible sequences Ca, 
(bj = 0, 1, N) . A feasilple sequence Ca is compatible 

with Mc, where there exists at least one triplet of 
corresponding numbers x, y, z, compatible with Mx, My, M^, 
respectively. ^ Since the list in Table 3 is exhaustive, it 

2S is possible to construct every feasible sequence Cm, and 
these sequences are used to search for solving triplets 
X, y, z, such that no valid solution is excluded from 
being found. 

30 
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List of ii"Values 



Recalling Restrictions (17)^ (18), and (19), it is 
seen that If, for some index me {0, W-1}, Cjb - 2 

then it is necessary that Q+i « 2 also* And if Qh-i =^ 2, 
then it is necessary that Cm ^ 0 (that is, C^i e {1, 2}). 
This implies that one of the following is true: 

1. There exists a boundary index n e {1, . . . , W) 
such that Cffl = 2 for all m > Ca ^ 1 for in « n - 

1 and Cb, e {0, 1} for all m < 

2. All of the carries are 2 {n = 0 in this case) . 

3. All of the carries are in {0, 1) (n = iff + 1 in 
this case) . 

A feasible boundary index n e {0, 1, W+1} is 

generally not unique, and there might exist several 
possible values for n. Therefore, construct a list of 
boundary index n-values which includes all of the values 
of n that correspond to solutions, and no other values of 
n. Clearly, for all n ^ m < N, = Q^+i =^ 2. So, looking 
in Table 3 it is seen that for all such m the pair 
(MCN,C») must be one of: 10,2), (1,2), (2,1), {3,1), 
(4,1), (5,1), (^,0), (7,0). Since C^-i « 1, Cx, » 2 for 
n 6 {1, W) it is inferred for such n that it is 

necessary for zn = n - 1 that the pair (MCN,asi) is one of: 
(0,0), (1,^2), (2,^1), (3,x), (4,^1), [5,x), [6,k) , 
(7, X), where ^1 means CN 1, 562 means OH ^ 2 and x means 
that CN may have any value. Additional restrictions on n 
are imposed by the mask Mc. It is necessary that Cn-i, Car 
are all compatible with this mask. 
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Given the masks, this permits the construction of a 
preliminary list of possible values of n. As seen below, 
this list is often too large, and some terms must be 
erased. 

Feasible Carry Seqfuence 

The sequence Cm is completed, given a value for n, by 
setting values to Co, Ci, These missing values 

of carries must all be in {0,1}. Hence, starting from 
Table 3 the list of pairs (Q,, Cxa+i) is modified by erasing 
from it all of the pairs which include 2. The remaining 
list, which is relevant to the construction of the 
missing carries, may be replaced by the following 
equivalent list of inference rules: 

• MCai = 0, CN = 1 + j' - Jt 
CN ^ -1: Ca = 0,C^i - 1 
CN = 0: = Q+i = 0 

CN « 1: Cja ^ Cja+i ^ 1 
CBI = 2: Cm ^ IfCjnn = 0 

• MCN - 1, CW = i + J 
CN = 0: Ga = 0 

CN = 1 : Cia = Cjni+i 

CN = 2: Ca = 1 

• MOT =2, CN = i - Jc or MCW ^ 4, CN = j - 

k 

cai = --1: Ca = 0 

CN = 0: Cj„ « Ca+i 
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• MCN = J, CN « i or MCN ^5, CS - j 

CN = 0: CjB = 1 => Cj^i == 1 

CN =^ 1; Cjff « 0 » CjB+i = 0 

CN=-1: = 1 => Cffl+i « 1 
CjM+i « 0 => C„ « 0 

CN« 0: C/a 0 « 0 • 

Qa+i = 1 =^ Cffi = 1 

• MCN « 7, CN « 0 

CN « 0: Nq restrictions 
Table 4. Carry Sequence Inference Rules 

like the list of pairs in Table 3 from which Table 4 
is derived, this set of inference rules is exhaustive in 
the sense that each feasible sequence Co, Ca-x, of 

{0/1} terms, must be compatible with these rules, and 
each such sequence, which is compatible with these rules 
and with Mcr is feasible. 

In setting values to the carries Co, Cn-i, these 

values are constrained by the mask Mc and by the inference 
rules of Table 4. In addition, C^-x 1 if 
n € {1, 2, M) , The mask Mc uniquely defines those 

terms of Ch which correspond to non ^x' -characters (note, 
however, that a ^2' character in Mc is permitted only for 
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m > n. Otherwise n should be erased from the list of 
jn-values) . The set of inference rules of Table 4 may be 
divided into three (not disjoint) groups: 

1. Assignment rules: for example^ Cju = Gn+i = 0; Gn = 
1; etc. 

2. Right continuation rules: for example, 
Q, = 1 Cffl+i = 1; Cffi = Cft+i; etc. 

3. Left continuation rules: for example, Cio+i = 0 
On ^ 0/ ^/B ~ Qj]+i; etc* 

In principle, the solution-seeking method 
illustrated in Figure 4 is used here also, where input 
list 401 contains a list of prospective n-values and 
block 411 operates as described below to search for a 
suitable carry sequence or output that there is no 
solution. 

First applying only the assignmen.t rules, it is 
possible to assign values to some of the carry terms. 
Note that there are several ways to deduce a definite 
value for a Qa (mask, assignment rules, Cn-i = 1) . It may 
happen that there may be contradictions. Therefore, each 
time a definite value is deduced for a given Qn, it is 
necessary to check to see if the particular carry bit was 
assigned a different value earlier. A contradiction means 
that this particular n should be erased from the list of 
n- values . 

Suppose all of the methods, described above for 
deducing a definite value for a C^, were used and no 
contradiction was found. Some of the defined carries may 
be neighbors (Cj^, Cj^+i). For each such pair of neighbors it 
is necessary to find the MCN and CN corresponding to the 
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index m, and test for a contradiction by the 
corresponding inference rule from Table 4. If there is a 
contradiction/ that zi must be erased from the list of 
n-values. If all of the pairs of neighbors were tested 
and no contradiction was founds the continuation rules 
are applied, one at a time. This process will create 
chains of consecutive defined carries*, separated by 
chains of consecutive (yet) undefined carries* As the 
process continues the chains of undefined carries shrink 
and it may happen that one of them disappears completely. 
That is, the right end of one chain of defined carries 
becomes a neighbor of the left end of the following chain 
of defined carries. Such neighbors must be tested for 
contradiction by the inference rules of Table 4, If any 
contradiction is found then that n should be erased from 
the list of /?~values- If the process ends and cannot be 
continued any further and no contradiction was found, 
then either all of the carries are defined and there is a 
complete, feasible, sequence of carries, or some chains 
of undefined carries were left over* In this case, a 
point was reached where no more contradictions are 
expected. It is then possible to choose one end of an 
undefined carries chain and choose for a value for that 
end of either 0 or 1, at random. No contradiction can 
arise from this operation, because, as was mentioned 
above, the set of inference rules of Table 4 is 
exhaustive. The new carry becomes a left or a right end 
of a chain of defined carries. The continuation rules are 
applied to this new end, again and again, until the end 
of the chain meets an end of another chain or until no 
further continuation rule can be applied, and then an 
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undefined carry is assigned at random. This process is 
repeated until all of the carries are assigned definite 
values. Note that if the new end meets another chain of 
defined carries, namely if the new end becomes a neighbor 
of another end there cannot arise a contradiction because 
the other end could not be continued at an earlier stage 
and this jaeans that its new neighbor may have the value 1 
or the value 0 without causing any contradiction. 

As was mentioned above, if any contradiction was 
found then the value of n must be erased from the list of 
73-values. If there are contradictions for all values of 
n, that is, if at the end, the list of n-values is empty 
then the fixed point generator must state that there 
exists no solution and stop. 

Through this process it is possible to discover if 
no feasible sequence of carries exists, and otherwise to 
produce, in principle/ every feasible sequence of 
carries. If there is a feasible sequence of carries it 
can be used to construct every triplet of solving numbers 
Yi Zi as described below. 

Number Construction from a Carry Sequence 

Note that at this point, all contradictions in the 
carry sequence have been eliminated. 

To construct the numbers x, y, and z, given a 
feasible carry sequence Cj^, first assume that this entire 
carry sequence is known and that it is feasible- For each 
value of in e {0, 1, ;^-l} there exist numerical 
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values for C^^ Cjn+i# MCN/ CN and perhaps some of the values 
1, jr k. 

Next, start with Equation (21) i + j - k 2Cja - Cjn+i 
and transfer to the right hand side of this equation all 
5 of the known values for i, j, k. This results in an 
equation of form 



10 



a = RBS 



(29) 



where the right hand side (RHS) has a known 
numerical value that equals a simple additive and/or 

subtractive combination (oD of the unknown values of i, 
15 j, k. 

It is easy to see that 
BBS = 2Gn - - CN * (30) 



20 and that a depends on MCN such that 

a + CN=i + j- Jc: (31) 
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For instance, if MCai » 3 then CN=i, a = j- ^ and 
Equation (29) becomes j - k ^ BHS. If MCN - 2 then 
CN = i - J^, a = j and Equation (29) becomes j = SHS, and 
so on. The form of Equation (29) and all of its 
5 solutions, in all of the possible cases, are suinmarized 
as follows: 

• MCKT = 0, a is an empty expression 
SHS =0: i, j, k are all known 

• MCN =1, a = -k 

10 RHS « -1 or 0: k ^ -RHS 

• MCN =2, a = J 

RHS= 0 or 1: j = RHS 



15 
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RHS - 1: a, j) - (0, 1), (1, 0) 
BBS = 2: (i, j) = (1, 1) 

RHS = -1: {1, j, k) = (0, 0, 1) 

KHS = 0: (2, j, k) - {0, 0, 0), (0, 1, 
1), (1, 0, 1) 

BBS - 1: (i, j, k) = (0, 1, 0), {1, 0, 
0), (1, 1, 1) 

BBS « 2: 11, j, k) - (1, 1, 0) 

Xable 5. Values for and k 

Thus/ knowing the nmn^rical values of MCU and RBS 
for every me. {0, 1, it is possible to select 

from Table 5 a solution which completes the triplet i, jf 
kg for every m. Wherever the list includes several 
solutions for some coiabination of MCN, RHS, one of the 
solutions is chosen at random. Making such choices for 
all values in e {0, 1, W-1} completes the 

construction of 

Fixed-point Generator Solution-Seeker 

Based on the above discussion, the solution-seeking 
method illustrated in Figure 4 can be used for the fixed 
point generator. First, list 401 is a list of possible 
n-values. Block 411 then operates as follows: In 
accordance with the procedures detailed above, try to 
construct the missing terms of C^. If there are any 
contradictions, there is no solution for the selected i2, 
and decision-point 413 branches to erase block 417. If 
there are no contradictions/ then sequence Cn is a 

IL 9-2001-0023 58 



feasible solution. Using the constructed sequence Cjb and 
masks 2»4, My, Mz, set values for i, j, and k for each 
m G {0, 1/ Whenever there is more than one 

possibility for choosing j, k, make a random choice. 
The found solution is the completed construction with x, 
Yf Zf and Co,, and decision-point 413 branches to output 
block 415 to return the solution. 
Operation of the Addition Slgnificand Generator 

Addition significand generator 235^ (Figure 2) always 
applies fixed point generator 245 with = p. The 
reasoning is as follows: the esiponent of ' a + ^ is almost 
always equal to the exponent of c. The only exception to 
this is in the post-normalization case which occurs when 
a + h rounds upward to produce - 1.00 — 0 exactly. In 
this case, the exponent of c is 1 greater than that of 
a + Temporarily ignoring this exceptional case (which 
is discussed below) , align the significands S^f S^t and Sc 
according to the values of and Qh- When this is done, 
some of the trailing bits of Sa, S/, are positioned to the 
right of the least significant bit of Sc and form ^^tails", 
as shown below, which contribute to the sum c only 
through carry and/or rounding into the least significant 
bit of the sum: 

Fixed point masks tails 

Sa: 0 

5i>: 0000 

Sci 

With reference to Figure 3, input 303 to fixed point 
generator 245 contains the Shifted and truncated masks of 
S^, as shown above. In addition, input 305 is a carries 
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mask Mc = ''Oxx. . .xCp", where Cp e {0, 1, 2} has the 
contribution of the tails, which combines the effects of 
carry and of rounding - 

First/ a numerical value is chosen for Cp of 0, 1/ or 
2, Next, the tails are generated, and finally Se and the 
left parts of and Sy, are generated. This is done by 
using the fixed point generator, as detailed below. 

Denote the leftmost bits of the tails by azr 
respectively, and the remainders of the tails by as, -bs* 

fixed point tails 
masks 



Sa: 0 [32 33 

Sjfc,: 0000 [b2 b% 

Scz 



10 Clearly, the value of Cp is determined by a^f b2, 33/ 

and Jba. 

Either Qa ^ {0^ 1} or Qi, e {0, 1}, and therefore 
either as ^ 0 or bs = 0. So, denoting cs = 23 + it is 
then clear that C3 = 33 if' fib e {0, 1} and C3 = bs if 
15 Da € {0, 1}. Cp thus depends on aa, bz, and 03, 



The Tails Triplet 



Actually it is not necessary to know the whole 
'20 sequence of the bits of C3, but only the result of an OR 
operation over all of the bits of C3, denoted herein as 
OR(c3)- The triplet of bits (32, ba/ ORic^s)) is herein 
denoted as the tails triplet". Thus, for instance, if 
the tails triplet is (1^ 1, 0) and the rounding mode is 
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round up then q, = 1. If the tails triplet is (1, 0, 1) 
and the rounding mode is round to nearest/even then 
Cp = 1, and so on. 

S All of the possible tails triplets, ' for each of the 

rounding modes and for each of the possible values of Cp, 
are as follows: 

• round down 

Cp = 0: (0, 0, 0), (0, 0, 1), (0, 1, 0), (0, 1, I), 

0 (1/ 0, 0), (1, 0, 1) • 

Cp = 1: a, 1, 0), (1, 1, 1) 

• round up 

Cp = 0: (0, 0, 0) 

Cp = 1: (0, 0, 1)\ {0, 1, 0)% (0, 1, 1)^ (1, 0, 

1 0)°, (1, 0, 1)^ (1, 1, 0) 

Cp = 2: (1^ 1^ 1)° 

• round to nearest/even 

Cp = 0: (0, 0, 0), (0^ 0, 1), (0, 1^ 0)o, (1, 0, 0)o 

Cp = 1: (0, 1, 0)°o, (0, 1, 1)^ (1^ 0^ 0)°o, (1, 0, 
1)S (1^ 1, 0), (1, 1^,1) 

Table 6. Tails Triple-ts 

Some of the triplets listed in Table 6 have a 
nximerical subscript and/or superscript: The subscript o 
means that the corresponding triplet implies round to 
even case. Such a case is possible with the indicated Cp 
value only If the last character of Se is forced to be 
*0'. The existence of a superscript indicates that the 
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rounding component of the ccmtribution of the tails is l, 
which means that a result with Sc = 1.00- is 
post-normalized and is potentially wrong/ because the 
exponent may have changed size/ thereby invalidating the 
original shifting assumption by 1 bit (as further 
discussed below) . 

The generation of S^, St, Sc starts by constructing a 
three-character ma^k for the tails triplet. The elements 
of this mask/ which correspond to az and jba are simply 
copied from Ms^t, Msb or are set to be '0^ if ag and/or Jbs 
fall outside of the range of the corresponding shifted 
mask. The element which corresponds to cs is set to be ^1' 
if the corresponding part of Msa or Msh includes at least 
one ^1' character. If not, the element will be ^0' if no 
^x' exists in the appropriate part of Msa or Msb and ^x' 
otherwise . 

After the mask of the tails triplet is ready, one 
tails triplet compatible with this mask is chosen from 
the complete list in Table 6, Note that in the case of 
round to nearast/even each of the triplets (0, 1, 0) / 
(1, 0^ 0) appears twice: once with Cp = 0 and once with 
Cp = 1- These appearances are considered to be distinct. 
Namely/ a choice such as (0, 1/ 0) with Cp = 0 is 
different from the choice (0, 1, 0) with Cp = 1. 

Now that a tails triplet has been chosen/ the 
construction of the two tails is straightforward- The 
construction of Sc and of the left hand parts of S^/ St is 
performed by the fixed-point generator. 
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If the tails triplet chosen from Table 6 does not 
have any superscript it means that even if the generated 

Sc is 1.00 0^ there is no post-normalization. In this 

case, the result is correct and acceptable. 

If^ however^ the tails triplet has a superscript and 
the fixed-point generator produces 5c 1.0C«0 then it 
means that some thing may be wrong with the result and it 
is not certain that the produced S^, S^, and Sc satisfy 
Sc « roiznd(2*^^S5 + Z'^Si,) , as it should. This is discussed 
here in further detail . 

If the superscript is ° then the result is definitely 
wrong and cannot be corrected. In such a case, the result 
should be discarded and the significands construction 
should be repeated - 

If the superscript is ^, it means that the result is 
in fact correct, and should be accepted. 

If the superscript is ^, it means that the C3 part of 
the tails should be constructed with care: Note that 
there are only two tails triplets with a superscript ^/ in 
the list: They are (0/ 1^ 1} and (1, 0, 1) in the round 
to nearest/evBn mode, with Cp - !• Since the resulting Sc 
was 1.00...0 it means that the p + 1 first bits of the 
exact sum were 011...1 and that C3 should be concatenated 
to the right of this, in order to produce the complete 
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exact stmt. This means that in order for the result 
generated by the fixed point generator to be usable / the 
leftmost bit of cs should be 1. The rest of the bits of C3 
are not important. If the leftmost bit of 03 cannot be 
chosen to be 1 because of mask constraints, it means that 
the solution should be discarded and significand 
construction should be repeated. 

Note that no possible solution of the significands 
generator problem is excluded by the method of generation 
described above, not even those vrith So = 1.00-,,0 and 
exact sum OH...IC3. 

Solution-Seeking Procedure for the Addition Significand 
Generator 

To use the solution-seeker illustrated in Figure 4/ 
it is first necessary to construct list 401, and this is 
done by using Q^, Qh and Msaf Msbi to produce a 
three-character mask for the tails triplet, and then 
constructing a sublist of tails triplets compatible with 
this mask from the complete list in Table 6. List 401 
contains these tails triplets. Block 411 operates by 
invoking the fixed-point generator to construct Sc and the 
left hand parts of Sa, S^,. If the fixed-point generator 
states that there is no solution, then decision-point 413 
branches to erase block 417. Otherwise, use the chosen 
tails triplet to construct the tails, and thereby 
complete the construction of Sq and -Sjj, after which 
decision-point 413 branches to output block 415, which 
outputs Sa and the zero-padded left hand parts of 5^, St 
as the found solution. 
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Operation of the Subtraction Slgiiificand Generator 

Denote c ^ a - bf and a - round (c) . The rounding 

error is denoted by s = I c - c I . If c is rounded down, 
c + 8 = a - or ^ + (c + e) a. If c is rounded up, 
then £ s = a - h, or b c ^ (a + s) . in either case 
there is an exact identity which includes only one 
addition of nonnegative numbers. The numbers can be 
considered to be fixed-point numbers, and if there were 
masks for these three fixed-point numbers, it would be 
possible to use the fixed-point generator to generate 
them. This is in fact the ca-se, as is shown as follows: 

Note that the non-zero bits of g always lie to the 
right of the least significant bit of a. Also, 
a ^ c ^ a > c, so the non-zero bits of s lie to the right 
of the least significant bit of a as well. This means 
that the bits of c + s are composed of the bits of a and 
the bits of s written in sequence, one after the other, A 
similar point holds for a + This can be illustrated 
graphically: Assime, for instance/ that the rounding mode 
is round down, that j2b = 3 and that Cc « 1- Then the masks 
for a, b, and c + s may be chosen to be 

Si, 

000 

So 

0 XX 

000 



bz 

c + e: 
a 
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These masks are composed of the masks for Ss, St, Scf 
padding '0' characters and 'padding ^x' characters. Note 
that/ because z can be anything, the mask for e has ^x' 
characters; furthermore, because a is the largest of the 
5 three numbers a, b, and c, a determines the shift of the 
other two. Because the rounding mode for this example is 
round down, the rightmost three characters for the mask Ma 
are all ^0'- If the rounding mode were round up, the 
rightmost three characters for the mask . would be ^x' 
10 instead. 

The fixed'-point generator is used to generate the 
three numbers and then extract from their binary 
representations the bits of Sa, Sj^, 5cr. In this particular 
case for the fixed-point generator, N - p ^ 3, and the 

15 mask for the sequence of carries is *^Oxx- . -xxO". It turns 
out that the maximum needed value for N is of the order 
of 2p, because if Si is larger than p it is necessary to 
know only if b > 0, and the details of the bits of b are 
of no significance. A similar treatment may be used if 

20 the rounding mode is round up. In such a case, however, 
the bits of e must be added to those of a instead of c. 
Details for the round down Roundinqf Mode 
The identity that must be used in the round down 
case is b + (c + s) -a. In a similar manner to Table 1, 

25 the combinations of Qb, Qc which should be considered are 
listed below in Table 7. 
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Qb 


0 0 11 


0 0 


. . 0 


0 


1 1 


, - 1 


Qc 


0 10 1 


2 3 


p-1 




2 3 


■ . P 





1 


2 3 




p-1 


^p 


2 3 


P 


>P 


Qa 


>P 


0 0 




0 


0 


1 1 


1 


1 








• 






« 







Table 7: Qi), Qc Fairs fox: Stibtractlon 

The coiranon length of the masks that must be 
presented to the fixed-point generator for most of the 
5 combinations is JfJ p + maxiOt^Qc) - The mask for Sj, should 
be padded to the left by Qb ^0' characters (if Qb > 0) and 
to the right by Qc - Oh ^0' characters (if Qb < Gc) - The 
mask for Sc should be padded to the left by Qo ^0' 
characters {if Qc > 0) and by Qb - Qa ^x' characters to 
10 the right (if Qc < ^) - The mask for Sa never has to be 
padded to the left (because a ^ b, ^ > c) and should be 
padded to the right by maK{Qb/Qc) '0' characters (unless 
Qi3 == Qc = 0) - The mask of the carries must always be of 
the form ^^Oxx. * .xxO". 

15 The cases {Qj^^Qc) 6 { (0/>p) , (l/>p) , (2p,0) , {>p,l) 1 

should be treated in a slightly different way: 

In the two cases where Qc ^ P ox Qc > Pt Sc - 0, 
unless this is not compatible with its mask, in which 
case there is no solution. Generating a solution, then, 

20 is straightforward. 
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In the two cases where Qb t p or Qh > Pt 5b may be 
any bit string which is compatible with its mask. The 
three-number masks that must be presented to the 
fixed-point generator are of length N - p + + 1 each, 
5 where the mask for h is composed only of ^0' characters 
except for the rightmost one, which is ^1' if Jb ^ 0 and 
is ^0' otherwise. The mask for So is padded by Qc '0' 
characters to the left and by a single 'x' to the right . 
The mask for is padded b^r + 1 ^0' characters to the 

10 right. After a solution to the fixed-point problem is 
determined, the bit strings for 5a/ St, and Sc may be 
easily constructed. 

Details for the round up Rounding Mode 

The identity that must be used in the round up case 

15 is b + c = + s) . The combinations of Q^, Qp which 
should be considered are also those listed in Table 7. 
Again, in most cases, the length of the masks, presented 
to the fixed-point generator is = p + jnax{Qi,,Qc) - The 
mask for is padded by Qh ^0' characters to the left (if 

20 Qb > 0) and by Qc - Qb ^0' characters to the right (if 
Qb < Qc) - The mask for Se should be padded by Qc '^O' 
characters to the left {if Qc > 0) and by Qb - Qc '0' 
characters to the right (if Qc < Qb) . The mask for Sa does 
not have to be padded on the left side. On the right side 

25 the mask for must be padded with Qc ^0' characters (if 
Qc > 0) and by Qb - Qc characters on the right (if 

Qc < Qb) • The mask for the carries is, again, of the form 

In the four cases where either Qc ^ p or Qb ^ p the 
30 treatment is, again, slightly different: 
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In the two cases where Qc ^ Sc must be zero, 
unless this is not compatible with the mask, in which 
case there exists no solution. 

In the two cases where Qe> > p, s = ^ and c = a 
(Oc = 1 is impossible, then) . So Si> may be chosen to be 
any number which is compatible with its mask and Sa - Sq 
may be chosen to be any positive number which is 
compatible with the masks of both Sa and of Sc. 

In the discussion above, it has been implicitly 
assumed that the exponents of c and of a are the same. 
This is always so in the case of round down. However, if 
roundi) is round up there exists one exceptional case: If 
the significand of c is 1.00,..0 and 8 > 0 then this 
implicit assumption is not satisfied. Unless the leftmost 
bit of g (the bit which corresponds to the leftmost ^x' 
character in the right padding of 5a) is 0, this leads to 
an error. So solutions returned by the fixed-point 
generator in which Sc = 1,00.,.0 and the leftmost bit of e 
is 1 should be rejected, and an additional attempt to 
produce a solution should be made. 

Details for the round to nearest/even Roundincr Mode 
For round to nearest/even the rounding is sometimes 
round up and sometimes round down. So the algorithm for 
this case is a mixture of the algoritlims for round up and 
for round down. Consider again the combinations of Qb^ Qg 
listed in Table 7. In the same way that the masks for S^, 
Si,, Sc were extended and padded for round down and round 
up, those masks are also extended and padded in the round 
to nearest/even case. A new factor, however, is the 
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splitting of the discussion of each combination of Q^^, Qc 

into four subcases: 

1- Case of round to nearest /down , in which the 
identity b + (c + s) - a is usedr and the extended 
masks of S^, St, Sc are padded in the same way as 
in the round down case/ except that the ^""kx^.^k'' 
padding that corresponds to s is replaced by 
^^Oxx padding* 

2. Case of round to nearest/up: The identity 
i + c « (a + e) is used, and the extended masks of 
S^, Sj>, Sc are padded in the same way as in the 
round up case, except that the padding, 
that corresponds to s is replaced^ again, by 
^'^Oxx...x" padding. As in the round up case a 
solution with Sc l.DO — 0 is rejected, unless the 
leftmost X (of s) is replaced in the solution by 0 
(i.e. the bits of e are compatible with *^00xx. ..x") 
or the bits of s in the solution are 0100. -.0 
exactly, 

3. Case of round to evBn/down: This is exactly like 
case 1 above, except that the ^*Oxx. ..x" padding is 
replaced by "100... 0'' padding, and the last 
character of the mask of Sq is replaced by ^0' 
(this can be done only if the original last 
character of Msc is ^0' or ^x' ) . 

4. Case of round to even/up: This is exactly like 
case 2 above, except that the "Oxx-.-X'" padding is 
replaced by "100., '.0'' padding and the last 
character of the mask of Sc is replaced by a ^0' 
(again, this can be done only if the original last 

IL 9-2001-0023 70 



character of Mso is '0' or 'x' ).'a solution with 
Sc="1.00...0 must be rejected. Thus, a solution 
with Se = 1.00..-0 cannot be produced in this case. 
Note, however, that such a solution may be 
5 produced in case 2 (see the discussion at the end 

of case 2, where the bits of b are 0100... 0). 
This discussion may be completed in a 
straightforward manner to also include the cases where 

^ p or ^ P- 
10 Solution-Seeking Procedure for the Subtraction 
Siqniflcand Generator 

Subtraction significand generator 240 (Figure 2) 
does not utilize sets of "tails triplets" as are employed 

I 

by addition significand generator 235, and therefore does 
15 not recursively use the solution-seeker method 
illustrated in Figure 4- Instead, the solution-searching 
action of subtraction significand generator 240 is part 
of the floating-point subtraction test-generator' s 
solution-seeking block 411 is as follows: Given numerical 
20 values for fib, Qc from calculator 225 based on the 
particular values of qt, qc selected from the 
floating-point subtraction test-generator's list 401, use 
Table 7 to classify this pair in order to use the 
appropriate procedure, out of those described above, to 
25 generate three significands. If the construction was 
successful decision-point 413 branches to output block 
415/ which outputs the complete floating-point solution 
(biased exponent and significand) • Otherwise, if the 
construction failed/ decision-point 413 branches to erase 
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block 417 which erases the selected qti qc selected from 
the floating-point subtraction test-generator's list 401. 
The Biased Exponent Generators 

Biased exponent generator 210 (Figure 2) includes 
5 both definite biased ejiponent generator 215 and 
indefinite biased exponent generator 220. In both of 
these biased exponent generators gi 6 {0/ 1} and 
©3 - ei + gi* There are then two possible cases: ea = ei 
and es - ei + 1. In the case of definite biased exponent 

10 generator 215, there is a given definite value of q2 such 
that ©3 = ^2 + ga ei + 5*1- That is, the exponent of the 
output result is a definite amount different from the 
exponent of either input operand. In the case of 
indefinite biased exponent generator 220, however, the 

15 value of qz is indefinite, and it is merely known that 
02 > g. That is, the exponent of the output result is not 
a definite amount different from the exponent of either 
input operand. The definite biased exponent generator and 
the indefinite biased exponent generator are non-limiting 

20 examples of a more general definite exponent generator 
and a more general indefinite exponent generator, 
respectively^ that obtain solutions for exponent pairs in 
these two cases. 

Definite Biased Exponent Generator 

25 In the case ©3 = ei, the common value of ©i and es 

must be compatible with both the masks Mi and M3* If the 
two masks have different number characters in the same 
position then they are incompatible and no pair ei, es 
exists. Otherwise, it is very easy to produce their 

30 intersection, M13, based on the individual character 
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intersections as shown in Table 8 below. Note that 
incompatible character intersections are not defined and 
are denoted by 0, 





'0' 










0 






0 








^0' 







Table 8. Mask Cliaractar Zntersactloixs 

The problem is now one of producing 82 and 03 that 
satisfy 63 = ea + gz^ There are masks for ez, ez (i.e. Mz, 
Miz) and it is possible to 'construct a mask composed of 
numerical characters only for gs- This is equivalent to a 
problem of fixed-point addition, and therefore, the 
problem may be solved by using the fixed-point generator - 

In the case 03 = ei + 1, note that the right hand end 
of the string of bits of biased exponent ei, must be one 
of the following: 0, 01^ Oil, 0111, ...,.Oll_.l (the last 
string is of length pit) . Because €3 = ei -f- 1, the right 
hand end of es must be, respectively: 1, 10, 100/ 1000, 
.../100,.^0 (here also the last string is of length t/) . 
Comparing the possible right ends of ei, es with the masks 
Mx, j% it is usually possible to erase some of the 
possibilities and what is left is a reduced list of pairs 
of right-hand ends of • ei, ea (which constitute 
solution-seeker list 401 in Figure 4) . In any of these 
pairs, the left ends of ei, ©3 must be identical. This 
means that the masks of ei and es may be chosen to be 
composed of known numerical characters in the right ends. 
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and of the intersection of the left hand ends of Mi and M3 
(if the left ends of Mi and Hz are incompatible then the 
corresponding pair of right ends, will be erased from 
list 401) . 

Thus, every choice of a pair of right ends of ei, ©3 
results in a condition similar to the one for the case ei 
- ©3: There are masks Jfe, M13 for es/ ©3 and a mask for qz, 
and ©2, ©3 must be found from the relation ©3 = e2 + qz. 
This, again, can be solved by the fixed-point generator. 
If the generator states that there is no solution it 
means that the selected pair of right ends will be erased 
from list 401, and another pair should be selected and 
tried. If list 401 is empty it means that there is no 
solution which satisfies es = es + 1. 

When employing the solution-seeker illustrated in 
Figure 4, given Hi, M3, and qz as a selection, the 
procedure for block 411 is as follows! 

1. If ii^ and M3 are incompatible, output that there 
is no solution. Block 411 is complete, after 
which decision-point 413 branches to block 417. 
Otherwise, if Mi and are compatible, construct 
M.3 and continue. 

2. If A&3 exists, present itfa, Mi3, and to the fixed 
point generator to generate ea and &3 satisfying 
©3 = ©2 + g2. 

3 If the fixed point generator states that there is 
no solution, output that there is no solution. 
Block 411 is complete, after which decision-point 
413 branches to block 417. Otherwise, if there 
exist ©2 and ©3 satisfying 63 » ©a + j^, construct 
©1 from ©3. 
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4. Return ei, 82^ and es. Block 411 is complete^ 
after which decision-point 413 branches to block 
415 to output the solution ei, ejy and as. 
Indefinite Biased Essponent Generator 
5 The analysis of the indefinite biased exponent 

generator is similar to the above analysis of the 
definite biased exponent generator, up to the point where 
there exists a new mask, W13, for ©3 (this applies to the 
case Bs = ei as well as to the case es ©a + 1) - Thus, the 
10 remaining problem is to generate q2f ezt and ©3 where 
there are masks Mz and for ea and ©3/ respectively, 
which satisfy the relation (63-62) = qfe > ^• 

The smallest B2 that is compatible with Mz is 
obtained by replacing each ^x' in Afe by ^^O', and this is 
15 denoted by e2 smallest • The largest 03 which is compatible 
with i^i3 is obtained by replacing each ""x^ in M12 by a , 
and this is denoted by as largest- There exists a solution 
for the indefinite biased exponent generator if and only 
if ©3 LAjtGEsi' - Bz stm^s^ > cj. If this inequality is not 
20 satisfied, the indefinite biased exponent generator 
reports that no solution exists ^ and is finished. 

If a solution exists, it is necessary to choose a 
random pair ez , ez for which e% - ez > q, and for which 
ez and es are compatible with 1^ and Mx3f respectively. 
25 This is done as described below; 

To start, erase from Mz all of the and '1^ 

characters, to leave a submask composed of ^x' characters 
only- The numbers ez'^ which are compatible with this 
submask, are in a natural isomorphic (one-to-one) 
30 correspondence with the numbers ez that are compatible 
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with Mz. Clearly/ ©2 is a monotonically-increasing 
function of ©2'' and vice versa- Similar relations exist 
between ea and via the mask Mx^. 

The construction of random ez and which are 

5 compatible with the masks Mz and M13, respectively/ and 
satisfy es' - ez' > q, is illustrated graphically in 
Figure 5 and is computed analytically as follows (in this 
description/ ez"", ©2** smallest/ ^2° lamest/ ^2*" correspond to 
©2/ ©2 SMALLEST/ ©2 LARGEST/ respect Ively / by the isomorphic 

10 correspondence via Mz7 and e^^, ez^ lAfJCssT/ ©3°' smallest 

correspond to 63/ ea largest/ ©3'/ ^3' smallest/ respectively/ 
via itfis) . Figure 5 shows an e2 axis 501 and an €3 axis 503. 
1. Compute Bz smallest/ ©3 largest with 62*" smallest - 00- . .0/ 
and ea* LARGEST = 11- --1/ respectively- A horizontal 
IS 02 smallest line 505 and a vertical B3 largest line 507 

thus delineate two boundary lines of the solution 
space/ which intersect at a point 509- For 
convenience/ also compute ez maximaii/ ®3 mimimal with 
^2** MAXIMAL = 11-. -1/ and 63° minimal = 00.. -0, 

20 respectively. A horizontal Bz maximal line 506 and a 

vertical 33 namt^ line 508/ and a point 510 at the 
intersection of line 506 and line 507 are useful 
in performing searches, as will be discussed 
below. 

25 2, Compute ez^ largest - maxiez^ I es laksbst - es > q*} . 

This may be done by any convenient 
search-and-test method, but it is well-known in 
the art that a convenient and rapidly-converging 
search merhod for monotonic functions is the 
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binary, or ^^bisection" search, where an interval 
having 28 elements is divided into two contiguous 

intervals of 6 elements each, whose boundary 
point can be quickly evaluated to determine in 
5 which of the two intervals the next iteration of 

the search should be performed (where an interval 
has an odd number of elements, 28 + 1, the 
'^bisection" breaks the interval into intervals of 
5 and 8+1 elements, respectively) . Graphically, 

10 this is represented on the line segment defined 

by point 509 and point 510, which may be 
successively bisected until ^2 iargiest is found. 
Also graphically, an eg largest point 511 
establishes an upper limit on the extent of the. 

15 solution space, A line 513 represents - ez - qz 

where 52 is the minimum value such that q2 > g*- 
Line 505;. line 507, and line 513 and the area 
enclosed thereby, thus define a solution space 
514, where ©3 - ^2 > g. Mote that it is not 

20 necessary to determine line 513 or the analytical 

counterpart thereof, nor is it necessary to 
determine all the solutions in solution space 
514- It is only necessary to determine a single 
point at random within space 514, such that every 

25 point within space 514 has approximately the same 

non-zero probability of being selected. 
3. Choose an integer e2*^* e [02^ sj^h^sTr j^^Ronsr] r at 
random, and from this derive ©2'- Graphically, 62' 
is represented by a horizontal line 518 through a 

30 point 515, which is chosen randomly in the line 
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segment defined by point 509 and point 511, 
inclusive of point 509 and point 511. 

4. Compute ©3°' suMiBST = min{ S3° | es - ea' > g}, and 
from this derive ea'.sMjaiissT. This can be done, for 
example, using a binary search, which is 
graphically represented on the line segment 
defined by point 515 and a point 516, which is 
defined by the intersection of line 508 and line 
518. A point 517 represents 63' sm&llesii' and it can 
be seen that point 517 necessarily lies on line 
513. 

5. Choose an integer 63*' e [ej,"' smcojssrr ©3* largest] at 
random, and from this derive 63'. Graphically, a 
point 519 chosen randomly on the line segment 
between point 515 and point 517 will have 
coordinates 62', ©3', which is guaranteed to be in 
solution space 514. 

Alternatively, it is possible to solve for ez"' and 
ea*', and afterward derive ^2' and 63'. Set the returned 
solution ©2 = 62' and 63 = ©3'. The right end of ©i is known 
and the left end may be copied from 63. Also, g2 ©3' - 
ez'. 

When employing the solution-seeker Illustrated in 
Figure 4, given Mi and M3 as a selection, the procedure 
for block 411 is as follows: 

1. If A?L and M3 are incompatible/ output that there 
is no solution. Block 411 is complete, after 
which decision-point 413 branches to block 417. 
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otherwise/ if Mi and Ms are compatible, construct 
Mi3 and continue. 

2. Using the above procedure, produce random eg' and 
es' which are compatible with and Wis/ 
respectively, and which satisfy es' - ©2' > <?• If 
such 62' and 63' do not exist, output that there is 
no solution. Block 411 is complete, after which 
decision-point 413 branches to block 417. 
Otherwise, if ez' and 63' exist, continue. 

3. Set ez « ez', ©3 = es', and = 63' - ea' and 
construct ei from ea. 

4. Return ©a, ©2, ^3, and gz- Block 411 is complete, 
after which decision-point 413 branches to block 
415 to output the solution ei, ez, es, and gjj. 

While the invention has been described with respect 
to a limited number of embodiments, it will be 
appreciated that many variations, modifications and other 
applications of the invention may be made. 
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